home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / dsycon.z / dsycon
Text File  |  1998-10-30  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSYYYYCCCCOOOONNNN((((3333FFFF))))                                                          DDDDSSSSYYYYCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSYCON - estimate the reciprocal of the condition number (in the 1-norm)
  10.      of a real symmetric matrix A using the factorization A = U*D*U**T or A =
  11.      L*D*L**T computed by DSYTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DSYCON( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK, IWORK, INFO
  15.                         )
  16.  
  17.          CHARACTER      UPLO
  18.  
  19.          INTEGER        INFO, LDA, N
  20.  
  21.          DOUBLE         PRECISION ANORM, RCOND
  22.  
  23.          INTEGER        IPIV( * ), IWORK( * )
  24.  
  25.          DOUBLE         PRECISION A( LDA, * ), WORK( * )
  26.  
  27. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  28.      DSYCON estimates the reciprocal of the condition number (in the 1-norm)
  29.      of a real symmetric matrix A using the factorization A = U*D*U**T or A =
  30.      L*D*L**T computed by DSYTRF.
  31.  
  32.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  33.      condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      UPLO    (input) CHARACTER*1
  38.              Specifies whether the details of the factorization are stored as
  39.              an upper or lower triangular matrix.  = 'U':  Upper triangular,
  40.              form is A = U*D*U**T;
  41.              = 'L':  Lower triangular, form is A = L*D*L**T.
  42.  
  43.      N       (input) INTEGER
  44.              The order of the matrix A.  N >= 0.
  45.  
  46.      A       (input) DOUBLE PRECISION array, dimension (LDA,N)
  47.              The block diagonal matrix D and the multipliers used to obtain
  48.              the factor U or L as computed by DSYTRF.
  49.  
  50.      LDA     (input) INTEGER
  51.              The leading dimension of the array A.  LDA >= max(1,N).
  52.  
  53.      IPIV    (input) INTEGER array, dimension (N)
  54.              Details of the interchanges and the block structure of D as
  55.              determined by DSYTRF.
  56.  
  57.      ANORM   (input) DOUBLE PRECISION
  58.              The 1-norm of the original matrix A.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDSSSSYYYYCCCCOOOONNNN((((3333FFFF))))                                                          DDDDSSSSYYYYCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      RCOND   (output) DOUBLE PRECISION
  75.              The reciprocal of the condition number of the matrix A, computed
  76.              as RCOND = 1/(ANORM * AINVNM), where AINVNM is an estimate of the
  77.              1-norm of inv(A) computed in this routine.
  78.  
  79.      WORK    (workspace) DOUBLE PRECISION array, dimension (2*N)
  80.  
  81.      IWORK    (workspace) INTEGER array, dimension (N)
  82.  
  83.      INFO    (output) INTEGER
  84.              = 0:  successful exit
  85.              < 0:  if INFO = -i, the i-th argument had an illegal value
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.